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I. INTRODUCTION 



During the last few years, the data collected from the acoustic peaks in the cosmic 
microwave background radiation [1] has allowed to obtain a more precise measurement of 
the baryon asymmetry of the universe (BAU). This is expected to further improve in the 
near future with the MAP experiment [2] and the PLANCK satellite [3]. At the present 
time, the measurement of the baryon-to-entropy ratio is 

0.7 x 1(T 10 < ^ < 1.0 x 1(T 10 , (1) 

where s = 2n 2 g if T 3 /45 is the entropy density and is the effective number of relativistic 
degrees of freedom. 

It is well known that in order to obtain an asymmetry starting from a symmetric state with 
a vanishing baryon number, three requirements must be satisfied: baryon number violation, 
C and CP violation, and departure from thermal and kinetic equilibrium [4]. In the standard 
model (SM) of electroweak interactions, the main source of CP violation comes from the 
phase Sqkm i n the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. Although 
this phase is able to account for the experimentally observed CP violation in the neutral 
i^-mesons and, as recently observed, in the Bd system, it has been shown that it is not 
possible to generate sufficient BAU through 5ckm [5]. Furthermore, the strength of the 
phase transition is too weak in the SM and the universe is approximately in equilibrium [6] 
(for reviews on electroweak baryogenesis, see for instance Refs. [7, 8]). 

In the context of supersymmetric (SUSY) extensions of the SM, it has been pointed out 
that the above problems can be in principle overcome [9]. Moreover, in the presence of 
light stops the electroweak phase transition can be strong enough for baryogenesis to take 
place [10-16]. Moreover, SUSY models contain new CP- violating sources beyond 5ckm, 
namely the Higgs bilinear term, /i, and the soft breaking terms (gaugino and squark soft 
masses, bilinear and trilinear couplings). These can be classified as flavor-blind or flavor- 
dependent. The first category includes the phases of the /i- and 5-parameters, of the 
gaugino masses and the overall phase of the trilinear couplings Aij. Two of these phases 
can be eliminated by U(1)r and U(\)pq transformations. The second category contains the 
phases of (after the overall phase is factored out), as well as the ones appearing in the 
off-diagonal elements of the soft squark masses. The low energy implications of these flavor- 
dependent phases on K and B meson CP- violating observables, as well as in rare decays, 
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have been extensively studied in the literature [17-19]. The question that naturally arises is 
whether the new SUSY phases could significantly enhance the CP- violating sources, so that 
super symmetric electroweak baryogenesis could account for the observed baryon asymmetry 
of the universe. 

A considerable amount of work has been done concerning the implications of the new 
flavor- independent CP- violating phases on generating an acceptable value of the baryon 
asymmetry [20-25]. It has been shown that in the minimal supersymmetric extension of the 
SM (MSSM), if the relative phase M between the gaugino soft mass and the \x term is not 
too small, M > 0.04, a considerable BAU can be generated through the scattering of the 
charginos with the bubble wall. However, the non-observation of the electric dipole moment 
(EDM) of the electron, neutron and mercury atom [26, 27] imposes severe constraints on 
the flavor-diagonal phases [28, 29], forcing them to be small. Since in this limit the theory 
does not acquire any new symmetry, one has to deal with a naturalness problem, and this 
is precisely the so-called SUSY CP problem. 

In particular, the EDM's bound the phase 0^ to be < 10~ 2 , if the SUSY particle masses 
are not too heavy (< 1 TeV). It has recently been claimed that new contributions to the 
EDM of the electron could eventually rule out the electroweak baryogenesis scenario based 
on flavor-independent CP- violating phases [30]. A possible way to generate enough BAU 
while evading these constraints is to work in the heavy squark limit {juq ~ 3 TeV) [31]. 

However, if we assume that SUSY CP violation has a flavor character such as in the 
SM [32], one is led to a scenario where all flavor conserving parameters as the /z-term and 
gaugino masses are real. In this framework, the dominant sources of CP violation that 
are relevant to electroweak baryogenesis are in general associated with the lightest of the 
right-handed up-squarks [33]. In the usual MSSM scenario with complex jj, and gaugino 
masses, this contribution would be always subdominant, and could even be neglected when 
compared to that of the charginos and neutralinos [21, 22]. 

In this paper we study the effects of the flavor-dependent CP- violating phases on the 
mechanism of electroweak baryogenesis. We show that in generic SUSY models with non- 
universal soft SUSY breaking terms, and in particular with non-universal A-terms, the 
squark contributions to the BAU are far from being negligible. Although in this framework 
the /i-term is real, flavor mixing might lead to a potentially large baryon asymmetry. As 
an example, we analyze the baryon asymmetry in supersymmetric models with Hermitian 
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flavor structures, a type of model that also provides an elegant solution for the EDM's 
suppression [32, 34, 35]. 

Obtaining sufficiently accurate transport equations for particles propagating in the pres- 
ence of a CP- violating bubble wall at the electroweak phase transition is crucial for the 
computation of the baryon asymmetry generated in the context of electroweak baryogenesis. 
In spite of the general consensus on the existence and nature of the CP-violating sources 
responsible for the baryon production, there is still a controversy in the literature in what 
concerns the strength and form of these sources as well as the transport equations to be 
used in the calculations [20-25]. For instance, the results obtained in the MSSM framework 
for the chargino and squark sources by making use of the continuity equations and the re- 
laxation time approximation [21, 22] are different from those based on WKB-methods [23]. 
Furthermore, it has been claimed in Ref. [23, 24] that the dominant contribution to the 
chargino source, typically found within the MSSM and which is of the form e^Vid^Vj (v i j2 
are the expectation values of the two Higgs doublets) [21, 22], is absent, thus leading to a 
suppressed baryon asymmetry [23, 24]. However, in a recent work [25] it is argued that such 
a suppression is in fact an artifact of the approximation used by the authors of Refs. [23, 24] 
in order to compute the CP-violating currents. In view of the above discussion, we adopt 
in the present work the approach of Refs. [20-22, 25] in order to derive the CP-violating 
currents and sources. 

The paper is organized as follows. By using the closed time path formalism, we build in 
Section II the quantum Boltzmann equations to obtain the CP-violating sources relevant 
to baryogenesis and discuss their general flavor dependence. In Section III we compute the 
baryon asymmetry induced by these CP- violating sources. An analytical approximate ex- 
pression for the asymmetry is also presented. In Section IV we obtain model independent 
bounds for the up-squark left-right mass insertions by requiring the baryon-to-entropy ratio 
to lie in the experimental range. We also consider some specific SUSY models with mini- 
mal and non-minimal flavor structure and discuss whether or not they can account for the 
required BAU. Our numerical results are presented in Section V. Finally, we present our 
concluding remarks in Section VI. 
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II. CP- VIOLATING SOURCES FOR THE BARYON ASYMMETRY 



Before computing the CP- violating sources for the baryon asymmetry, we shall briefly 
recall some characteristics of non-equilibrium quantum field theory [36, 37]. During a first- 
order phase transition, the thermodynamical system is far from equilibrium. In order to keep 
an explicit time dependence, one should use the real-time finite temperature quantum field 
theory. The most used and powerful formalism to describe such a system is the so-called 
closed time path formalism (CTP), which is a generalization of the time contour integration 
to a closed time path. More precisely, the time integration is deformed to run from — oo 
to +00 and back to —00. The main effect of this closed time path is to double the field 
variables so that for each field we have four different real-time propagators on the contour. 
In the case of a boson field 0, we can write the corresponding Green functions in terms of 
the 2x2 matrices [38] 

G*=\ * t I , (2) 



G> -G t 



where 



G>{x,y) = -i{<t>{x)<l>\y)) , G<(x, y) = -i (0%)#r)> , (3) 
Gl(x,y) = 9(x,y) G>(x,y) + 6(y,x) G<(x,y) , 

G%x,y) = 6(y,x) G>(x,y) + 9(x,y) G<(x,y) . (4) 

In what follows the subscript will be omitted to simplify our notation. 

As mentioned in the introduction, an accurate computation of the CP- violating sources 
responsible for electroweak baryogenesis is crucial to obtain a reliable estimate of the baryon 
asymmetry. Here we adopt the method developed by Kadanoff and Baym [36] to derive the 
quantum Boltzmann equations for a generic bosonic particle asymmetry. We shall compute 
the sources using the formalism based on CTP and the Dyson equations as described in 
Ref. [21]. 

In general it is possible to write the Dyson equation for any propagator G as 
G(x,y) = G°(x,y) + J d A Xl J d A x 2 G°{x, Xl ) E(x 1 ,x 2 ) G(x 2 ,y) 

= G°(x,y) + J d A Xl J d A x 2 G(x,xi) t(x u x 2 ) G°(x 2 ,y) , (5) 
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where G° are the non-interacting Green functions and S are the self-energy functions, which 
can be also expressed in terms of the 2x2 matrices 



£> — E* 



satisfying the relations 



j/) = 0(a;, y) S > (x, y) + 6>(y, x) S < (x, y) , 
£*>,y) = 0(y,z) S>(x,y) + 0(z,y) S<(x,y) . (7) 

Iterating Eq. (5) one finds 

G(x, y) = G°(x, y) + J d 4 x x J d 4 x 2 G°(x, Xl ) t{x u x 2 ) G°(x 2 , y) 

+ J d 4 xi J d 4 x 2 J d 4 X 3 J d 4 x^ G°(x, xi) £(xi, x 2 ) G°(x 2 , rr 3 ) S(x 3 , x±) G°( X 4,y) 
+ .... (8) 

In order to compute the sources for the squark diffusion equations, we shall write the 
quantum Boltzmann equations (QBE) for a generic bosonic particle asymmetry, 

where is the number density of particles minus antiparticles and S is the associated 
CP-violating source. More specifically, we shall derive the QBE for the current 

(Jf(x)) = i (4>\x) <f>(x)) = (n+(x),J*(x)) , (10) 

or in terms of the Green functions, 

(■/£(*)) = - lim(^ - G<(x,y) . (11) 

x >y 

It proves convenient to work in the coordinate system of the center of mass, which is 
defined by 



Hence 



X = (T,R) = ^(x + y), x = (t,r)=x-y. 



G<(x,y) = G<(t,r,T,R) = -i(tf(T ~\,R-\) ^T+ t -,R + r -)) . 



The next step is to find a solution for G < when the system is not in equilibrium. This can 
be achieved by applying the operator ( □ +m 2 ) on both sides of the equivalent representations 
of the Schwinger- Dyson equation (5): 

(D x + m 2 ) G(x, y) = 5 4 (x, y)l + J d 4 Xl t(x, Xl ) G(x 1: y) , (12) 
G(x, y) (D y + m 2 ) = 5 4 (x, y) 1 + J d 4 Xl G(x, x x ) t(x u y) , (13) 

where 1 is the identity matrix. Since (SG)< = S* G < -S < G* and (GE)< = G* E<-G< £*, 
the equations for the G K component read as 

(n x + m 2 )G < (x,y) = J d 4 x, [^(x, Xl ) G < (x 1 ,y) -E<(x, Xl ) G\x u y)] , (14) 
G<(x,y) (n y + m 2 ) = J d 4 x 1 [G t (x,x 1 )E < (x 1 ,y)-G < (x,x 1 )J:\x 1 ,y)]. (15) 

The variation of the current is thus given by 
d J$(X) 



OX" =-^{^-^)G<(x,y) 

= - J d 4 x x [E*(X,xi) G < {x 1 ,X) -S<(X,x 1 ) G\x u X) 

- G t (X,x 1 ) E < (x 1 ,X) + G < (X,x 1 ) E { (x u X)] . (16) 

Using Eqs. (4) and (7), we can compute the CP- violating source for a generic squark 
current: 

T 

S= 9 U ^ + V ■ J^X) = - J d 3 r' J dt' {T>(X,x') G < (x',X)-Z < (X,x')G > (x',X) 

— oo 

- G > (X,x') ^<(x',X) + G < (X,x') S > (x',X)} . (17) 

Inserting the interactions of the <fi field with the background, the self-energy functions 
S > (x, y) can be written as follows 

E>(x,y) = g(x) 5\x - y) + g(x) G 0> (x,y) g(y) + / d 4 z \g(x) G%x, z) g(z) G°(z,y) g(y) 



+ Jd A w J d 4 z [g(x) G°(x, w) g(w) G°(w, z) g{z) G°(z, y) g(y)] > + ... , (18) 

where the scalar function g(x) describes the interaction of the field at position x with the 
background fields. 



Taking into account the chiral and flavor structures of the squarks, this generic squark 
source can be written as 

T 

Sf = ~\ dV J dt> { [g(X)]f 5\X - x>) 5f [G<(x>, X)]f 

-[9(X)]f 5\X-x')5? B [G>(x',X)]f 

-[G>(X,x')]f [g(x% B 8\X-x')5f 

+ [G<(X,x')\t [g { x% B 5\X-x')5? B 

+ [9(X)]f k C [G>(X,x')\c [g(x% B [G<(x',X)]f 

-[9(X)]f [G<(X,x')f k [g(x% B [G>(x',X)]f 

-[G>(X,x')]f [g(x')]f [G<(x>,X)f k [g(X)f k f 

+ [G<(X,x')]f [g(x')]f [G>(x',X)f k [g(X)f k B + 0([g(x)f)} , (19) 

where (i, j, k) are flavor indices and (A, B, C) refer to L, R chiralities. In the above expres- 
sion, we have kept only terms up to second order in g(x). Note that this formula is valid for 
squarks of any flavor and/or chirality. Thus, we have 

< J£ = Sab ^ Sf B = , (20) 

with 0f = Uii and 0f = Um, the left and right-handed up-squark fields, respectively. As 
expected, the contributions of the terms proportional to <5 4 (X — x') vanish 1 . The structure 
of the source can be symbolically depicted as in Fig. 1, where the ® denotes interactions 
with the background fields, parametrized by the scalar function g(x). 

Before proceeding, let us consider the Higgs-squark interaction terms, arising from both 
SUSY conserving and SUSY breaking terms, that can be parametrized in the following way 
(flavor indices are omitted) : 

C 9 = Ul [g(X)} LL U L + U* R [g(X)] RR U R + (U* L [g(X)] LR U R + H.c) . (21) 

Note that in addition to the interactions with the background, the functions [g(X)] AB contain 
flavor and chirality mixing terms that originate from the soft breaking SUSY Lagrangian and 
become relevant when one computes higher order terms like those appearing in Fig. 1(b). 

1 Throughout we are using the relations: G < {x) = G > (—x) , [G < (x)]* = —G < (—x) , [G > (x)]* = 
-G>(-x) , which imply C7 < (a;) = [G>(x)]*. 
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(a) (b) 
FIG. 1: A symbolic representation of the structure of the CP-violating squark source. 



Within a generic supersymmetric extension of the SM, the functions g(X) can be 
parametrized as follows: 



[g{x)t- = (m| + vl(x) k h T u + Dl(x) 
[g(x)\T = (mI r + v\{x) hi k + dUx) 



'j 



[g(X)]ff = {- Vl (X) vK + v 2 (X) (Y U A )*).. , 

[g(X)] RL = [g{X)] LR] . (22) 

In the above equations, h u denotes the Yukawa coupling matrix for the up quarks and v 1,2(2^) 
are the x-dependent vacuum expectation values of the MSSM Higgs fields H\ 2 , defined as 
Vi — v cos P/y/2 and t>2 = wsin/5/v^2, with tan/3 = v 2 jv\. Mq and My R are the squark soft 
breaking mass matrices, /i is the Higgs bilinear term and Y A the soft trilinear terms, which 
are decomposed as 

Qtf) y = {h u ) l3 Al , (23) 

with no summation over 

The D-term contributions are given by 



Dt(X) = ^- - - sin 2 6 w j {v f {X ) - V *(X)) , 

Dl{X) = Q sin 2 9 w ^j 9 -±f- [v\{X) - v\{X)) , (24) 
g, g' are the SU(2) and U(l) weak couplings, respectively. These terms can be neglected 
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at the electroweak phase transition since they are suppressed by the weak couplings when 
compared to the terms proportional to the top Yukawa coupling. 

Since we are interested in performing a basis independent calculation, we should express 
the current as an invariant quantity. Let us define the invariant current as 

"/ 

d i J ^ = E sAA = Tr sAA > ( 25 ) 



i=i 



where n/ is the number of up-squark flavors. 

Under a generic rotation that transforms the left- and right-handed up squarks as Ul 
W L U L , U R ^> Z r U r , we have 



[G(x,y)>><] R ^ Z R [G(x, !/)>'<] fl 4= \G(x,y) 



>>,<i l 



(26) 



[g(x)] LL - W L [g{x)] LL Wt = [g(x)] LL 
[g(x)] LR - W L [g(x)] LR 4 [g(x)] LR 
[g(x)] RR - Z R [g(x)] RR Z R ee [g(x)} 



RR 



(27) 



Thus the invariant source which is given by 

± n J: 



Tr S AA = - f dV f dt'J2 {l9(X)ti C [G>(X,J)]° [g(x')f kt A [G<{x\X)\f 

- lg(X)]f k C [G<(X,x')f k [g(x')] C kl A [G>(x',X)\t + (X <- x>) + . . .} , (28) 
can be rewritten in any weak basis. The above equation can be recast in the form 

T 

Tr S AA = -2i j d 3 r' I dt' { Im aO]£ [G < (x , ,X)]f) 

oo *> fe=1 

x {[9{X)\f Wt? - W)\tk l9(x)]Z A ) +•••}• (29) 

Finally we obtain 

T 

Tr S RR = 4 I d 3 r' I dt' ( Vl (X) v 2 (x') - v 2 (X) v 1 (x')) { Im {[ G< ( x ^')\k [G < (x',X)]?) 

i,k=l 

— oo ' 

x Im (jm [Z R (Y A f Wl] [W L K 4] J } . (30) 
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A convenient basis for the current associated with the right-handed up-squark Um is that 
where the Green functions are diagonal. For simplicity, we shall assume that the left-handed 
squark mass matrix M LL is already diagonal. Therefore one has Wl = 1. Moreover, we 
will also assume that the left-handed squarks are heavy {tuq ^> m t ) and nearly degenerate, 
so that Mll — diag(mg, toq, tuq) and G(x)% = G(x) L for any flavor. The right-handed 
squark mass matrix M RR is diagonalized by a unitary matrix Z R so that 

Z R M RR Z R = d 2 RR , (31) 

d RR = &\&g{m R1 ,m R2 ,m m ). In this case, 

T 

nf 

TiS RR = A j d 3 r' I dt'(v 1 (X)v 2 (x')-v 2 (X)v 1 (x'))J2{lm([G < (X,x')] L [G < (x',X)]f) 

x lm^[z R (Y u A f h* u Z R ]J] . (32) 

In order to obtain an analytical expression, it is useful to perform an expansion in the 
bubble wall velocity v w . Such an expansion is well justified in the case of the MSSM, 
since bubbles are typically formed with thick walls (L w ~ (10 — 100)/T) and propagate with 
extremely nonrelativistic velocities (v w ~ 0.1 — 0.01) [39, 40]. Using the derivative expansion 

^) = E^(^T^)(^-^T, (33) 

it is easy to check that the first nonzero contribution to Eq. (32) is given by the n = 1 
term in the expansion (33). The associated contribution to the source is proportional to the 
function 

Vl (X) d» v 2 (X) - v 2 (X) d» Vl (X) = l - v\X) d» x (3{X) , (34) 

which in turn implies that the source (32) is linear in the wall velocity v w . Neglecting terms 
of higher order in v w , which amounts to using the thermal equilibrium Green functions, one 
can obtain an explicit expression for the dominant contribution: 

Tr S RR = v\X) d x (3{X) |E J ^ Im [ Z * K 4] J | , (35) 
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with I RR given by [21] 

oo 

Prr = J dk 4n*u Q u* [( 1 + 2Re (^)) IWr^qJq) 
o R 

+ (l + 2Re«)) 7(o;Q,r Q ,^,P i? ) + 2(lm«)+Im(n Q )) G(u r ,T r ,u q ,Tq)] . 

(36) 

The squark equilibrium distribution functions are 

n\ Q = 1/ [exp (u^q I T + iT\ Q / T) - l] , (37) 

where the finite widths T R q account for the interactions with the plasma, 

u 2 Q = k 2 + m 2 Q + m 2 Q (<P,T) , ^ = fc 2 + m^ + m^(0,T) , (38) 

and rhQ Ri ((p,T) are the field-dependent contributions to the squark masses, which include 
the temperature dependent self-energies HQ jRi (T) [41]. The functions I and G can be written 
as follows: 

I(a, b, c, d) = I+(a, b, c, d) + I-(a, b, c, d) , G(a, b, c, d) = G + (a, b, c, d) + G-(a, b, c, d) , 

(39) 

1 1 / a±c 

I±(a,b,c,d) = —-. — — — — — — sin I 2arctan 



2(a±c) 2 + (b + d) 2 \ b + d 

G ± (a, 6, c, d) = ^ (a±c)2 | (fe + rf)2 cos ( 2 arctan . (40) 

Before proceeding to the explicit computation of the baryon asymmetry, it is worth 
emphasizing that from Eqs. (35) and (36), the dominant contribution to Tr S RR will be 
associated with the lightest of the right-handed squarks, typically £r. In this framework, 
usually known as the light stop scenario, to satisfy the out-of-equilibrium conditions the 
mass of the right-handed stop should be small enough so that at the electroweak phase 
transition temperature T ew one has m^ R (T ew ) ~ m t (T cw ) = /i t f 2 (T ew ). This is precisely the 
case we shall consider from now on. In this context, the right-handed soft breaking masses 
are chosen so that m 2 R1 ~ m 2 R2 ~ mQ,m R3 = m 2 R uiq. In particular, the parameter m 2 R 
could be negative provided that no charge or color breaking minima appear [10]. 

In this case, 

Tr S RR ~ v 2 (X) d x p(X) I RR Im [z R (Y u A f h* u Z R ] J . (41) 
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In the limit where the trilinear terms are flavor conserving, one can easily recognize the 
well-known expressions obtained in Refs. [21, 22]. 



III. ELECTROWEAK BARYOGENESIS 

In the present scenario, the baryon asymmetry generated in the broken phase is deter- 
mined by the density of left-handed quarks 2 , created in front of the bubble wall in the 
symmetric phase [20-23]. Such densities induce weak sphalerons to produce a nonvanishing 
baryon number [42-44]. If the system is near thermal equilibrium and the particles are 
weakly interacting, the particle densities rii are given by rii = kifiiT 2 /6, where fa are the 
local chemical potentials and ki are statistical factors equal to 2 (1) for bosons (fermions) 
and exponentially suppressed for particles with masses m 8 ^> T. Assuming the supergauge 
interactions to be in thermal equilibrium and neglecting all Yukawa couplings except those 
corresponding to the top quark, it is possible to express til in terms of the densities of the 
chiral supermultiplet Qi = (q,q) , n,Q i = n qi + n qi , and of the right-handed top quark and 
squark, T = t R + t R . Indeed, assuming that all the quarks have nearly the same diffusion 
constant, one obtains due to the strong sphaleron processes, uq 1 = uq 2 = 2(uq + tit), 
implying the relation 

n-L = n Ql + n Q2 + n Qz = 5n Q + 4n T . (42) 

We can write down a set of coupled diffusion equations for the relevant densities uq , 
fix , tih = nu x + n H2 and = uh 2 — ■ Ignoring the curvature of the bubble wall, all 
the quantities become functions of z = r + v w t, the coordinate normal to the bubble wall 
surface. In the rest frame of the bubble wall we obtain 

D q n% - v w n' Q - TyN, - 6T SS N 2 - T m N 3 + lQ = 0, (43) 
D q r% - v w n' T + T y Ni + 3T SS N 2 + T m N 3 - lQ = 0, (44) 
D h n" H - v w n' H + TyN, - T h ^- = 0, (45) 

Kh 

D h nl - v w n' h + P r,iV 4 - (T h + 4r^)^ = 0, (46) 



2 In principle, the baryon asymmetry is determined by the density of all left-handed fermions, including 
leptons. However, in the usual electroweak baryogenesis scenario, there is essentially no lepton asymmetry. 
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where, according to Eq. (41), the CP- violating squark current is given by 

7Q ~ v w v\z) P'(z) I RR Im [p [Z R (Y u A f h* u Z R ] J . (47) 

Moreover, 

n Q n T n H + pn h n Q n T n Q + n T 

7Vi — , iv 2 - z- — h y : , 

kg k T k H kg k T k B 

AT U Q UT M n Q HT n « + n hlP (AA \ 
iV 3 — — , 7V 4 — , (48 J 

k Q k T k Q k T k H 

with the parameter p varying in the range from to 1. The coefficients k B = 3 , kg = 6 , 
kr = 9 , k H = 12 are statistical factors and the quantities 

T y ^ 2 ^thl a s T * 2.4 x 10~ 2 T , 

T ss = 16 a 4 s T ~ 3.3 x 1(T 3 T , 

T m = ^ sin 2 /3 ~ 4.8 x lO" 2 T , (49) 

r ft = — ^ ~ 7.1 x 10~ 3 T , 
140 T 

are reaction rates: T y corresponds to the SUSY trilinear scalar interaction involving the third 
generation squarks and the Higgs Hi plus all SUSY and soft breaking trilinear interactions 
arising from the superpotential term h t H 2 QT, T ss is the strong sphaleron rate, I\ and T rn 
arise from the Higgs and axial top number violating processes, while V ^ corresponds to the 
Higgs bilinear term. The numerical estimates of these rates are obtained assuming the top 
quark Yukawa coupling h t — 1, the strong coupling constant a s = 0.12, v(T) ~ T and 
tan (3 = 10. Finally, the diffusion coefficients are given by 

D q ~^, D h ~™. (50) 

Let us assume that T y and T ss are fast enough so that T y D h , T ss D q , T y D q » i> 2 . In 
this case Ni = OiYy 1 ), N 2 = 0(TJ^) and 

n Q n T n H + pn h _^ ^g n T | g ng + nr „ Q (51) 
kg k T k H kg k T k B 

This implies 

kg(9k T -k B ) ( v ( . 

k H (k B + 9k Q + 9k T ) 
k T (9k T + 2k B ) 

Ut = ~~T7l — Tlvi r~o/ — ( n n + P n h) , 53 

k B {k B + 9kg + 9k T ) 



14 



and 



n H + pn h l-p 2 n h 
N 3 ~ , N 4 ~ JVi — . (54) 



Thus, Eq. (46) can be rewritten as: 

D h nl - v w n' h + pTyNi - [(1 - p 2 )r, + T h + 4r M ] ^ = 0. (55) 

First we notice that for typical values of the scattering rate due to the top quark Yukawa 
coupling [20], the interaction rate (1 — p 2 )T y is fast enough in a wide range of values of p 
(except those values very close to 1), leading to the solution 

n h ~ 0. (56) 

Indeed, assuming v w ~ 0.1, the condition (1 — p 2 )T y D h > v 2 w together with Eqs. (49)-(50) 
imply that p < 0.998. 

Next, from Eqs. (52), (53) and (42) we find: 

k H (k B + 9k Q + 9k T) 

HH * k Q (9k T -k B ) ^ (57) 
k T 9k Q + 2k B 

•* - ~T Q HE^TtT "* (58) 

9k Q k T - 5k Q k B - 8k T k B 
= M9^ - k B ) (M) 

Taking a linear combination of Eqs. (43)- (45), which is independent of the fast rates T y 
and T ss we obtain an effective diffusion equation for uq : 

Dti'q - v w n' Q - Tn Q + 7 = 0, (60) 

where 

L> = - [D q {9k Q k T + k Q k B + Ak T k B ) + D h k H (k B + 9k Q + 9k T )] , (61) 

r = kB + 9k £ + 9kT (r m + r h) e( z) , (62) 

^ k Q{ 9k T -k B ) iqQ{z1 (63) 
A = 9k Q k T + k Q k B + Ak T k B + k H (k B + 9k Q + 9k T ) , (64) 

and 0(z) is the step function, which accounts for the fact that the rates T m and Th are active 
only in the broken phase. 
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To estimate the baryon asymmetry we need to solve Eq. (60) in the symmetric phase 
(z < 0). Imposing the boundary conditions n<g(±oo) = and the continuity of uq and ti'q 
at z = 0, a simple analytical approximation can be given: 



n Q (z) = A Q e^ D , 



where 



oo 



±^vl + AYD 
2D 



(65) 



(66) 



Assuming that the weak sphalerons are inactive in the broken phase (i.e. inside the 
bubble), n B will be constant in this phase. To find the value of this constant, one has to 
solve the diffusion equation for the baryon asymmetry in the symmetric phase: 



v w n' B = -0(-z) [n F T ws n L (z) + R n B (z)\ , 



(67) 



where rip = 3 is the number of families, R = \npY ws is the relaxation coefficient [45] and 
T^s ~ 120 a b w T is the weak sphaleron rate. The above equation is easily integrated and we 
find 

o 

-n F Y % 



n B 



(68) 



Using Eqs. (59) and (65), one gets 

9kQk B — 8k B kr — 5/cr/c, 



n B = —ufY v 



bkq 



DA 



Q 



(69) 



k Q (9k T - 9k B) \DR + vl 

To obtain a reliable approximation for the squark current, it is necessary to specify the 
Higgs profiles as functions of z, i.e. the functions v(z) and (3(z) which appear in Eq. (47). 
A simple expression is given by the kink approximation [22] 



v(z) = l -v{T) 



P(z) = P - -AP 



. . 2a 
tanh I a — —z 

-t-'ll! 



2a 

1 + tanh ( a — —z 



(70) 
(71) 



where L w /(2a) parametrizes the thickness of the bubble wall, a ~ 3/2 and A/3 is the 
variation of the angle fi(z) along the bubble wall, which lies in the range 10~ 2 > Aj3 > 
4 x 10 -3 for 100 GeV < m A < 200 GeV [39]. At first order in v w , all the dependence in z 
of the source 7q is given by the combination v 2 (z)f3'(z) . Substituting the profiles (70) and 
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(71) into Aq (cf. Eq. (65)), it is possible to integrate explicitly the z dependence. Finally 
we have 



An = 



v\T)Ap k Q (9k T -k B ) / r T f] x 



with I RR given by Eqs. (36)- (40), 

h + = 2a(a ~^ )( 6 a ~ 2) e"** 2 F!(1, 3 - |; 4 - |; -e- 2 ") + a(a - 4)(a - 2) tt e""* cosec(avr/2) 
- 2(1 + e 2a y 3 e~ 4a [(a 2 + 8)e 4a + 2a(a - l)e 2a + a(a - 2)] , (73) 

a = \+L w / (2a) and 2-^1(0, c; 2) is the hypergeometric function. 

In Section V we will present our numerical results for the baryon asymmetry using the 
approximate solution given by Eqs. (69), (72) and (73). 



IV. BARYON ASYMMETRY IN SUPERSYMMETRIC MODELS 

In this section we study the flavor-dependent contribution to the squark current, which 
enters in Eq. (41) through the imaginary part 



Tr S RR oc Im (jm [z r (Y u A f h* u Z R ] J . (74) 



33/ 

Our goal is to maximize the above quantity and, simultaneously, satisfy the experimental 
constraints on the off-diagonal A terms as well as on the CKM mixing matrix. 

To establish a direct connection between Srr and low energy supersymmetry phenomenol- 
ogy, let us express Srr in terms of the up left-right mass insertions, 5^ R . 

normalized to the mean value of the squark masses; U L,R are the unitary matrices 

that diagonalize h u , 

U L h* u U R ^ = d u = — diag (m u , m c , m t ) . (76) 

For the sake of simplicity and to maximize Tr Srr , let us assume that Z R — 1 . We notice 
that the unitarity of Z R implies that any deviation from the identity (or its permutations) 
will introduce additional suppression factors in the CP-violating source. From Eqs. (74) 
and (75) we obtain 

Tr Srr oc ^h 2 t \^L i m [ Ut fuR (^y . (77) 
17 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Im < S LR> 3 , 

FIG. 2: Baryon asymmetry as a function of Im (fi^^si for different values of m^. We assume 
/i = 500 GeV, v w = 0.04 and (m~) ~ uiq = 1 TeV 2 . The other parameters are indicated in Table I. 
The dashed lines correspond to the lower and upper bounds of the observed baryon asymmetry 
given in Eq. (1). 

Assuming for instance that the dominant contribution comes from large mixing between 
the t- and either the u- or c-quark, i.e. taking = U-^ = l/y/2 (i — 1,2), we have 

Tr S RR ~-^ t ^Im[(5 u LR y 3l ] . (78) 

In Fig. 2 we plot the baryon asymmetry as a function of Im (5^)3i for different values 
of rriA- It is seen that even for light pseudo-scalar masses m^, large values of Im {S^^)^ are 
required in order to reproduce the observed BAU. In Fig. 3 the baryon asymmetry is given 
as a function of Im (S^ R ) 3i for different values of fi. We notice that it is possible to get a 
sufficient BAU if 

> 0.14 , for /i^lTeV. (79) 

It is worth noticing that these bounds on the imaginary part of the mass insertions 
{S^ R )3i, 2=1,2 are compatible with the bounds obtained from the chargino contributions 
to Bd — Bd mixing and the CP asymmetry of Bd — > J/ipKs [46]. In addition to fulfilling 
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FIG. 3: Baryon asymmetry as a function of Im (#5^)3? for different values of [i. We assume 
vtla = 150 GeV, v w = 0.04 and (m~) ~ vtuq = 1 TeV 2 . The other parameters are taken as in 
Table I. The dashed lines correspond to the lower and upper bounds of the observed baryon 
asymmetry given in Eq. (1). 

the bounds on the LR mass insertions that arise from satisfying the BAU requirements (as 
shown in Figs. 2 and 3), the textures for the Yukawa couplings and trilinear terms must 
be compatible with the FCNC constraints [47], charge and color breaking bounds [48] and 
those associated with the measurements of the electric dipole moments of the neutron and 
mercury atom. Regarding the latter, the current experimental bounds are [26] 

d n < 6.3 x 10~ 26 e cm , d Hg < 2.1 x 10~ 28 e cm . (80) 

These values can be translated into bounds for the imaginary parts of the up and down 
(Slr) 11,22 [29]. In particular, 

Im (5*) £B < 10~ 6 and Im (5^) LR < 10~ 7 - lO" 8 , (81) 

from the neutron and mercury atom EDM's, respectively. 

Besides the underlying EDM problem that is associated with the imaginary parts of the 
diagonal elements of 5^ R , one should also bear in mind the extensive array of constraints 
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on the non-diagonal entries. These stem from FCNC bounds, limits on rare decays and the 
observed amount of CP- violation in the K and B d meson systems [49]. 

We can also express Eq. (74) in a different manner using the definition (23). Assuming 
\i real and keeping only terms proportional to h t , one has 

Tr S RR oc /A 2 \U 3 L k \ 2 Im [(Z R ) 3j A^U^U*^] , (82) 

where a summation over the indices j, k, I = 1, 2, 3 is understood. In this case, the imaginary 
part of Tr Srr is given by the simple form: 

Tr S RR oc /A 2 | U 3 L k 1 2 I U£ 1 2 Im (A u k3 ) , (83) 

where Zr — 1 has been assumed. The important point is that Srr is proportional to the 
imaginary part of A k3 , which are quantities weakly constrained by experimental data in 
contrast to other terms proportional to A^ = 1,2) [47]. 

In the case of hierarchical Yukawa couplings, the matrix U R can be assumed to be very 
close to the identity and the maximum values for U L are bounded by the CKM matrix. In 
such scenario, we have 

Tr Srr oc vh\ | U 3 L k 1 2 Im (A u k3 ) < rf* \ Vg KM \ 2 Im (A u k3 ) . (84) 

It is easy to see that the contribution proportional to A± 3 and A% 3 is very suppressed by \ V td \ 
and \ Vt. s \ 2 and thus, the dominant contribution comes from the imaginary part of A 33 which 
is always constrained in GUT-inspired SUSY models by the EDM bounds. In particular, in 
the minimal supergravity inspired model where the A-terms are universal and the overall 
phase is constrained to be < 10 _1 , this contribution becomes quite negligible and one cannot 
obtain a viable BAU through the squark current. 

Therefore, the flavor- dependent phases will play an important role in baryogenesis only 
if significant mixing between the top quark and up and/or charm quark is present. This 
mixing can be obtained in supersymmetric models with non-universal soft breaking terms 
and, for instance, with nearly democratic Yukawa couplings. Of course, this means that 
processes like t — > u + 7, + 7 may have much larger branching ratios than those 

predicted in the SM. 

To illustrate our results, let us now consider some simple textures for the quark and 
squark mass matrices at the electroweak scale. In case A we present the generic case of a 
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SUSY texture where the CKM matrix is dominated by the down quark mixing. Case B 
corresponds to SUSY models with Hermitian Yukawa couplings and trilinear terms at GUT 
scale. Of course, the RGE running down to the electroweak scale will induce some deviations 
from hermiticity, but their contributions will not be relevant to our BAU analysis. In case C, 
the textures are chosen such that the BAU produced at the electroweak scale is maximized 
and the phenomenological constraints are still satisfied. 

Case A 

We first consider the following simple texture for the trilinear matrix A u : 





fl 


1 


1 


\ 


A u = A, 


1 


1 


1 






V 


1 


a e i7r/2 


J 



with a = 0(1) a real parameter. We also assume that the up quark Yukawa coupling matrix 
h u is diagonal, i.e. h u = d u , Ul — 1 , Ur — 1. In this case one obtains 

Tr S RR oc fih 2 Im (A% 3 ) = aA tf ih 2 . (86) 



Case B 

Let us now assume that the up- and down-quark Yukawa couplings as well as the tri- 
linear ones are Hermitian matrices. In this case the up-quark Yukawa coupling matrix is 
diagonalized by a unitary matrix U such that h u = U T d u U* . Since the trilinear matrix A u 
is Hermitian too, then Im A 33 = 0. The squark source Srr can be expressed as 



Tr S 



RR 



■fih 2 t | U 33 1 2 { I U 31 1 2 Im ( A u 13 ) + | U 32 1 2 Im ( A u 23 



.)} • 



(87) 



It is easy to see that in order to maximize the imaginary part of Srr, a large mixing 
between the third and the first (or the second) families of up quarks is needed. For instance, 
assuming maximal mixing between the c and t quarks, i.e. 



U = 



\ 



1 
l/y/2 -l/y/2 
^0 l/y/2 l/y/2 ) 



(88) 
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and the following Hermitian texture for the trilinear matrix A u , 



A, = A< 



1 1 

1 1 b e in ' 2 
1 be"™' 2 1 



\ 



(89) 



(90) 



with b a real parameter of order 0(1), one obtains 

Tr S RR oc ^fxh 2 Im (A^) = ^A tf ih 2 . 

The corresponding texture for the down-quark Yukawa matrix h d is fixed by the structure 
of the CKM matrix, V C km = UD^ , where D is the mixing matrix that diagonalizes the 
down quark mass matrix, h d = D T d d D*, with d d = — diag (m d ,m s ,m b ). Thus, we have 



Vi 

h d = U T V£ KM d d V£ KM U* 



(91) 



Case C 

In this case we take 



1 1 ce"' 2 
1 1 1 
1 1 1 



U T 



( 



\ 



1 

1 

1 



where c = 0(1). The squark source Srr is then given by 

Tr S RR oc fxh 2 Im (A^ 3 ) = cA t fxh 2 



u R = t, 



(92) 



(93) 



Let us remark that in the previous examples most of the entries in the trilinear matrix 
A u were set for simplicity equal to one. It is clear that in general those elements can be 
nondegenerate and even complex. 

In the next section we shall present some numerical examples for the textures considered 
above and compute the BAU generated by the corresponding CP-violating squark current. 



V. NUMERICAL RESULTS 

In order to produce the required BAU, it is needed not only to violate the CP and C 
symmetries, but the universe has also to be out of equilibrium during the stage of baryogene- 
sis, to avoid any washout by the electroweak sphalerons. This can happen if the electroweak 
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phase transition is of first-order. To freeze the action of the sphalerons, their rate in the 
broken phase has to be smaller than the Hubble rate, i.e. 



r ew (T c ) < H{T C ) , 

where T c is the critical temperature of the phase transition. Using the sphaleron rate in the 
broken phase [50], this equation is translated into a condition on the vacuum expectation 
value of the Higgs fields at T c : 

^>l. (94) 

C 

In the standard model of electroweak interactions, this would imply a Higgs mass lighter 
than 43 GeV, which is already experimentally ruled out [27]. A way to avoid this constraint 
is to add scalar fields (with Higgs field-dependent masses) such that their contribution to 
the finite temperature potential will increase the strength of the phase transition. Naturally, 
the MSSM is an appealing alternative to the SM. In particular, it has been shown that the 
presence of a light right-handed top squark with small mixing considerably enhances the 
strength of the phase transition [10-12]. However, in this case a few constraints need to 
be satisfied. Within the MSSM, all these requirements impose restrictions to the allowed 
parameter space [16]: (i) a heavy pseudoscalar mass m A and a large tan/? regime, m A > 
150 GeV, tan (3 > 5; (ii) heavy left-handed stops with a relatively small mixing, rriQ > 1 TeV, 
0.25 < (A t - fi cot (3) jrriQ < 0.4; (iii) a light right-handed stop, 105 GeV < m~ tR < 165 GeV. 

In the context of models with non-universal A terms, we can expect to have more freedom 
and that some of these conditions may be relaxed. In fact, if we want to keep alive the light 
stop scenario, the constraints imposed by the electroweak phase transition (cf. Eq. (94)) 
are typically the same than in the one squark generation case, except that now we should 
require 

lfc)d< 0.4x^^ 0.07, (95) 

\ m ql 

if one assumes that (m~) = ttIq = 1 TeV 2 . 

From the last equation and the bound given in Eq. (79) (see also Figs. 2 and 3), we 
conclude that it is difficult to produce enough BAU using the flavor-dependent phases of 
the trilinear soft breaking terms unless > 1 TeV. 

In Fig. 4 we present the baryon-to-entropy ratio n B /s as a function of the \i parameter for 
different up-quark Yukawa coupling and trilinear term textures (cases A, B and C considered 
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TABLE I: Input parameters used in our numerical calculations [20-22]. 
MSSM PARAMETERS 

m w = 80.4 GeV m z = 91.2 GeV v = 246.22 GeV 

m u = 2.3 MeV m c = 0.68 GeV m t = 174.3 GeV 

A t = 400 GeV m Q = 1 TeV m R = 

m H = 115 GeV tan (3 = 10 100 < m A < 200 GeV 

BAU PARAMETERS 

T = 100 GeV c/* = 125.75 

L w = 20/T a = 3/2 4 x 10~ 3 < A/3 < 10~ 2 

Statistical factors 

kg = 3 kjj = 12 

kQ = 6 Izt = 9 

Diffusion coefficients 

A, = 6/T = 110/T 

Reaction rates 

r M = 0.1 t r ws = 5.4 x 10~ 6 T 

T ss = 3.3 x 10 -3 T T y = 2.4 x 10~ 2 T 
T m = 4.8 x 10~ 2 T T h = 7.1 x 10" 3 T 

Tq = Tr = Yj = 0.12 T 

in Section IV). The coefficients a, b, c that appear in the textures (85), (89) and (92), are 
chosen so that the CP- violating squark source is maximized and the constraints coming from 
the electroweak phase transition are satisfied. We also require the lightest Higgs mass to be 
consistent with the present experimental lower bound, m# > 109 GeV [27]. In particular, 
to maximize the CP- violating source and, simultaneously, satisfy the bounds on the lightest 
right-handed squark mass, mi R > 105 GeV, we take a = 1, 6=1 and c = 1.8 . From the 
figure it is seen that textures A and B cannot generate enough BAU. This is related to the 
fact that the coefficients a and b cannot be increased without violating the bounds on mi R . 
It is also clear that the produced BAU in the Hermitian case (B) is further suppressed by a 
factor 1/4 due to the Yukawa quark mixings. On the other hand, texture C can produce the 
observed baryon asymmetry provided that the \i parameter is large enough. We also notice 
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FIG. 4: The baryon-to-entropy ratio ub/s as a function of the \x parameter for different up-quark 
Yukawa coupling and trilinear matrix textures (cases A, B and C considered in Section IV). We 
have chosen a = 1, 6 = 1 and c = 1.8 for the coefficients in the textures (85), (89) and (92), 
respectively. We assume tjia = 150 GeV and v w = 0.04 . The rest of the parameters are chosen 
according to Table I. The dashed lines correspond to the lower and upper bounds of the observed 
baryon asymmetry. 

that in this case the coefficient c is less constrained by mi R . 

The baryon asymmetry as a function of the wall velocity v w is plotted in Fig. 5 for different 
values of fi and assuming the textures of case C, which succeeds in producing enough BAU. 
As we can see from the figure, to get a BAU compatible with the observational limits on 
n B /s (cf. Eq. (1)) large values of /i (/i > 700 GeV) are required. 

We should emphasize however that the electroweak baryogenesis is a strongly out-of- 
equilibrium process and all our computation is based on several approximations. Unfortu- 
nately, it is very difficult to estimate the errors done during the calculation because typically 
they come from different sources. Nevertheless, we can conclude that the scenario with a 
light stop and CP violation coming from flavor-dependent phases is in general disfavored. 
Of course, the constraints coming from the electroweak phase transition can be in principle 
relaxed in extensions of the MSSM with new scalar fields, for example by adding singlet 
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FIG. 5: The baryon-to-entropy ratio hb/s as a function of the wall velocity v w for different values 
of /i and rriA = 150 GeV. The other parameters are chosen as in Table I. We assume the textures 
given in case C of Section IV. The dashed lines correspond to the lower and upper bounds of the 
observed baryon asymmetry given in Eq. (1). 

fields (NMSSM) [51]. 

In the present model, where the CP- violating sources arise from the flavor- dependent 
phases of the SUSY soft breaking trilinear terms, the strongest constraint comes from the 
lightest up squark mass. Indeed, two conditions that will push the lightest up-squark mass to 
smaller values must be fulfilled. First, tur < T ew , otherwise the contribution of the lightest 
right-handed squark to the BAU will be exponentially suppressed by a Boltzmann factor. 
Secondly, as we have seen from Fig. 3, Im (5£ R ) 3i > 0.14, which typically means that some 
of the Ajj's have to be of the order or larger than rriQ. In other words, the 6x6 up-squark 
mass matrix will have large mixings. 

In this paper, we have restricted ourselves to the range of parameters satisfying the upper 
bound (95) such that all the experimental and theoretical constraints on the squarks and the 
lightest Higgs are satisfied in the light stop scenario. An important question that remains 
is whether it is possible to find GUT patterns for the trilinear terms, SUSY soft-breaking 
squark masses and quark Yukawa couplings, such as to obtain large values of 5^ R and a strong 
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first-order phase transition, while still satisfying the experimental limits on the squark and 
Higgs masses, as well as the EDM's constraints. 

VI. CONCLUSIONS 

Recent EDM bounds impose severe constraints on the usual scenario of supersymmetric 
electroweak baryogenesis based on flavor-conserving CP- violating phases. A natural solution 
to this problem is to work in the framework where all the flavor-conserving parameters, such 
as the //-term and gaugino masses are real. In this case, the dominant contribution to the 
baryon asymmetry is associated to the flavor-dependent CP- violating squark sources. 

In this work we have studied in detail the impact of non-universal A terms on the scenario 
of electroweak baryogenesis. By generalizing the standard approach, we have obtained the 
expression for the CP- violating squark sources with explicit flavor dependence. We have 
shown that if we impose on these terms the condition to have a strong first-order phase 
transition induced by a light right-handed squark, the baryon asymmetry of the universe 
produced at the electroweak scale is typically too small, thus disfavoring this scenario. 

On the other hand, if we assume that the problem of the strength of the first-order elec- 
troweak phase transition is solved through another mechanism as it can happen in extensions 
of the MSSM with additional Higgs scalars, it is possible to have textures for the up- and 
down-quark Yukawa coupling matrices in order to maximize the CP- violating source. This 
however implies a large mixing between the top quark and one of the light up-quarks (u or 
c) as well as large deviations from universality for the A terms (typically, (S^ R ) l3 23 > 0.14). 
Such large S1 R could have important implications on flavor-changing top decays. 

Supersymmetric electroweak baryogenesis is an attractive mechanism to explain the ob- 
served baryon asymmetry of the universe. Not only the physics involved in this process is 
directly related to low-energy observables, but it can also be testable in accelerator experi- 
ments in the near future. In particular, searches for a light Higgs boson and a light stop at 
LHC and Tevatron will constitute a test of the viability of this scenario. 

There are still a few questions to be answered and controversial issues to be clarified. For 
instance, the precise details of the electroweak phase transition are still unknown and there 
is at present a debate regarding the structure of the CP- violating currents that are relevant 
to baryogenesis. Nevertheless, a considerable progress has been done during the past few 
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years, and the effort directed to resolve these problems could give us definite answers in the 
near future. 
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